Collective oscillations of two colliding Bose-Einstein condensates 



P. MaddaloniB, M. Modugno, C. Fort, F. Minardi, and M. Inguscio 
INFM - European Laboratory for Non linear Spectroscopy (LENS) and Dipartimento di Fisica, Universita di Firenze, 

Largo E. Fermi 2, 1-50125 Firenze, Italy. 
(February 1, 2008) 



O ' 

o : 
o . 

(N ■ 

. 



(N 
> 

O 

m 
o 
o 
o 



i 

c 

o 
o 



Two 87 Rb condensates (F = 2,m/=2 and m/=l) are pro- 
duced in highly displaced harmonic traps and the collective 
dynamical behaviour is investigated. The mutual interaction 
between the two condensates is evidenced in the center-of- 
mass oscillations as a frequency shift of 6.4(3)%. Calculations 
based on a mean-field theory well describe the observed effects 
of periodical collisions both on the center-of-mass motion and 
on the shape oscillations. 
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Since the first realization of Bose-Einstein condensa- 
tion with dilute trapped gases 0], systems of conden- 
sates in different internal states have deserved attention 
as mixtures of quantum fluids. In this context, the im- 
portant issue of the interaction between two distinct con- 
densates was early addressed at JILA Q with the pro- 
duction of two 87 Rb condensates in the hyperfme levels 
\F = 2,nif = 2) and |1, — 1) in a Ioffe-type trap. Sub- 
sequent experiments of the JILA group have focused on 
the dynamics of two condensates in the states |2, 1) and 
1 1,-1) having nearly the same magnetic moment, con- 
fined by a time-orbiting potential (TOP) trap. In these 
experiments |^-||] the authors have investigated the ef- 
fects of the mutual interaction in a situation of almost 
complete spatial overlap of the two condensates. The re- 
sulting dynamics reveals a complex structure, and it is 
characterized by a strong damping of the relative mo- 
tion of the two condensates |3). More recently another 
group j^] has experimentally investigated a mixture of 
87 Rb condensates in different m/ states in a TOP trap, 
but no effects of the mutual interaction have been ob- 
served. 

These experiments have raised several interesting ques- 
tions from the theoretical point of view, such as the origin 
of the above mentioned damping [Q, the phase coher- 
ence properties of the two condensates || and the nature 
of Rabi oscillations in presence of an external coupling 
H , thus providing a challenge to explore new dynamical 
regimes in different experimental configurations. 

In this work we demonstrate an experimental method 
for a sensitive and precise investigation of the interac- 
tion between two condensates. These are made to col- 
lide after periods of spatially separated evolution and we 
get quantitative information from the resulting collective 
dynamics. Following the scheme originally introduced 
by W. Ketterle and co-workers [Eol to create an atom 



laser out-coupler, we use a radio-frequency (rf) pulse to 
produce two 87 Rb condensates in the states |2,2) = |2) 
and |2, 1) = |1). Due to the different magnetic moments 
and the effect of gravity, they are trapped in two poten- 
tials whose minima are displaced along the vertical y axis 
by a distance much larger that the initial size of each 
condensate. As a consequence the |1) condensate, ini- 
tially created in the equilibrium position of |2) , undergoes 
large center-of-mass oscillations, in a regime very differ- 
ent from that explored in ||] and analyzed in 0, where 
the two condensates sit in overlapping traps. The fact 
that the two condensates periodically collide opens the 
possibility of detecting even small interactions through 
changes in frequency and amplitude of the oscillations. 
Indeed, the periodic collisions of the |1) condensate with 
the 1 2), initially remaining almost at rest, strongly affect 
the collective excitations of both the condensates: (i) the 
center-of-mass oscillation frequency of the |1) condensate 
is shifted upwards; (ii) the shape oscillations of |2) con- 
densate, triggered by the sudden transfer to the |1) state, 
are significantly enhanced. 

The complex dynamics is quantitatively analyzed and 
found in agreement with the theoretical predictions de- 
rived by the numerical solution of two coupled Gross- 
Pitaevskii (GP) equations at zero temperature. 

We prepare a condensate of typically 1.5 x 10 5 87 Rb 
atoms in the F — 2, m/ — 2 hyperfme level (|2)), con- 
fined in a 4-coils Ioffc-Pritchard trap elongated along the 
z symmetry axis pd| . The axial and radial frequencies 
for the 1 2) state, measured after inducing center-of-mass 
oscillations of the condensate, are u> Z 2 — 2n x 12.6(2) Hz 
and lv±2 — 27r x 164.5(5) Hz respectively, with a mag- 
netic field minimum of 1.75 Gauss. By applying a rf 
pulse, the initial |2) condensate is put into a coherent 
superposition of different Zccman \rrif) sublevels of the 
F = 2 state, which then move apart: |2) and |1) are low- 
field seeking states and stay trapped, |0) is untrapped 
and falls freely under gravity, while | — 1) and | — 2) 
are high-field seeking states repelled from the trap. All 
the condensates in different Zeeman states are simultane- 
ously imaged by absorption with a 150 fis pulse of light 
resonant on the F — 2 — > F' — 3 transition, shone 30 ms 
after the switching-off of the trap. By fixing the dura- 
tion and varying the amplitude of the rf field, we con- 
trol the relative population in different Zeeman sublevels 
Q. A 10 cycles rf pulse at 1.24 MHz with an ampli- 
tude B r / = 10 mG quickly transfers ~ 13% of the atoms 
to the |1) state without populating the |0), | — 1) and 
| — 2) states. The |1) condensate experiences a trap- 
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ping potential with lower axial and radial frequencies 
(uji = U02/ v2) and displaced along the vertical y axis by 
yo = g/ujj_ 2 ~ 9.2 jum. After the rf-pulse, the |1) conden- 
sate moves apart from |2), and begins to oscillate around 
its equilibrium position. Due to the mutual repulsive in- 
teraction, the latter starts oscillating too, though with a 
much smaller amplitude. The resulting periodic superim- 
position modifies the effective potential, which is the sum 
of the external potential (magnetic and gravitational) 
and the mean-field one. We have studied the dynam- 
ics of the |1) condensate in presence ("interacting" case) 
and in absence ( "non- interacting" case) of |2), by varying 
the permanence time in the trap. We restrict to perma- 
nence times so short, i.e. less than 40 ms, that we can 
neglect both atom losses due to the condensate finite life- 
time, 0.7(1) s, and the heating (dT/dt = 0.11(2) fiK/s). 
For the non-interacting case, we have used a stronger rf- 
pulse (B r /=20 mG) that completely empties the |2) state 
coupling all the atoms in |1) and in the other untrapped 
Zeeman states, which rapidly leave the trap. 

We start considering the case of the |1) condensate 
in absence of the |2) condensate (non-interacting case). 
In Fig. [j]a we plot the center-of-mass oscillations as 
a function of the trapped evolution time in units of 
— 0.967 ms, after 30 ms of ballistic expansion. 
The center-of-mass undergoes sinusoidal oscillation with 
a measured frequency of cjxi = 27r x 116.4(3) Hz. This 
value is, as expected, a factor \[2 smaller than 
Furthermore the experimental data give no evidence of 
damping within the maximum monitored trapping time 
(26 ms). Fig. [If} shows the center-of-mass evolution 
of the two condensates in the interacting case. The 
center-of-mass of the |1) condensate exhibits a substan- 
tially different behaviour, which allows a clear quanti- 
tative analysis. The oscillation frequency is up-shifted 
to luj_i = 2tt x 123.9(3) Hz, i.e. a 6.4(3)% change with 
respect to the non-interacting case. This can be under- 
stood considering that the mean-field repulsion of the 
second condensate produces an anharmonic correction 
to the effective potential experienced by the |1) conden- 
sate. Furthermore, the oscillations appear now damped 
with an exponential decay time of about 60 ms. Indeed, 
each time the two condensates superimpose (about every 
8 ms), there is an energy transfer from the |1) conden- 
sate toward the |2) condensate. As a consequence we 
expect, and we do observe, effects on both the center-of- 
mass motion and the collective shape-oscillations of the 
1 2) condensate. 

A quantitative description of these experimental fea- 
tures requires the solution of two coupled Gross- 
Pitaevskii (GP) equations. Neglecting the interaction 
with the thermal cloud, the two condensates evolve ac- 
cording to 
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i = 1,2, where Vi are the trapping potentials: 

V 1 (x,y,z) = ^u;l 1 [(x 2 + y 2 )+X 2 z 2 ] (2) 
V 2 (x, y, z) = [(x 2 + (y- y f) + A 2 z 2 ] (3) 

and the asymmetry parameter is A = lo z /uj± ~ 0.0766 
for both traps. For the 87 Rb scattering lengths we use 
«22 = a\2 = 98 ciq and an = 94.8 ao 

Our experimental configuration allows to simplify 
these equations by using the fact that we are in the 
Thomas-Fermi (TF) regime due to the large number of 
atoms, and that the system is strongly elongated along 
the z axis, A <C 1 p4[ - For an elliptic trap, the low- 
frequency excitations with m = are linear superpo- 
sitions of the monopole (n = 1, I = 0, m = 0) and 
quadrupole (n = 0, I = 2, m = 0) modes p5| , ^6[ . The 
dispersion laws for the two modes, at leading order in A, 
are given by 
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In this limit the two frequencies are quite different, and 
the axial and radial collective excitations are almost de- 
coupled. The radial width is characterized by small- 
amplitude oscillations with frequency uji ow modulated by 
a large-amplitude oscillation with frequency Whigh, and 
vice versa for the axial width |Tij] . Therefore, since the 
interactions mostly affect the radial motion, we assume 
the axial dynamics to be still characterized by the low 
frequency oscillations of the TF regime. Then, we study 
the trapped dynamics in the x, y plane by solving the GP 
equation (Jl|) , by approximating our system as a uniform 
cylinder JOI. 

We start from an initial configuration corresponding 
to the stationary ground-state of Eq. (0), with all the N 
trapped atoms in the |2) condensate. Afterwards, at t = 
0, N\ atoms are instantaneously transferred from the \2) 
to the |1) state, the former remaining with N% = N — Ni 
atoms. Here we consider N% = 0.13A for the interacting 
case, and N\ — N for the non-interacting one. 

The theoretical curves in Fig. |l| show an up-shift of the 
center-of-mass oscillation frequency for the |1) conden- 
sate occurring in the interacting case, as experimentally 
observed. This shift is of 5.4%, in qualitative agreement 
with that measured. Furthermore, in presence of inter- 
actions the model correctly predicts a damping, which 
is not due to any dissipative process (the total energy is 
conserved), but to a transfer of energy from the center- 
of-mass oscillations of the |1) condensate to the other 
degrees of freedom of the system j2^] . Still, this damping 
time is nearly a factor 2 longer than that experimentally 
observed. 

To understand the origin of the discrepancies between 
theory and experiment, it's worth discussing the main 
approximations of our model. First, since the model is 
basically 2-dimensional, the energy transfer in the axial 
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direction is overlooked. Secondly, we completely disre- 
gard the interaction between the two condensates during 
the expansion p2| . Indeed, in our experiment, during 
the switching-off of the trap the two clouds acquire dif- 
ferent velocities and cross each other in the fall, but we 
expect that they are so dilute that we can neglect their 
mutual mean-field repulsion. Finally, our model doesn't 
take into account the elastic scattering, occurring when 
the relative velocity of the two condensates exceeds the 
sound velocity [p3[ . This effect, whose description lies 
beyond the mean-field approximation, represents an im- 
portant channel of both atoms and energy losses which 
plays a significant role, for example, in the four-wave 
mixing experiments with Bose condensates [ p4| , p5| . 

We consider now the |2) condensate. In the interacting 
case, both small center-off-mass oscillations and signifi- 
cant features for the aspect ratio oscillations emerge from 
our model. The latter are initially induced by the sudden 
change in the internal energy, consequent on the trans- 
fer of Ni atoms from the |2) to the |1) level In 
Fig. U we compare the theoretical evolution of the |2) 
condensate aspect ratio, i.e. the radial to axial width 
ratio, in absence and in presence of collisions. In the for- 
mer case, a faster (radial) oscillation superimposes to a 
slow (axial) one as a consequence of the decoupling be- 
tween the two oscillation modes. The frequencies uJhigh 
and ijJiow were separately measured by means of reso- 
nant modulation of the trapping magnetic field Jll[] . The 
non-interacting behaviour predicted in Fig. [2] should of 
course yield when there is only one trapped state. For 
instance this is the case of an atom laser out-coupled in a 
single-step transition, as for sodium |lO|j , or for rubidium 
in F = 1 state p7| . In the interacting case, we see from 
Fig. H that a significantly different behaviour is expected. 
Before the first collision at t ~ 8 ms, the oscillations are 
rather small since they are determined only by the sudden 
change in the number of atoms. Instead, at longer times, 
the changes in the aspect ratio become more pronounced 
due to the energy transfer during collisions between the 
two condensates. Once the ballistic expansion is taken 
into account, the simulation results are compared to the 
experimental data in Fig. ||. The agreement is only of 
a qualitative character. We attribute the discrepancy 
to the approximations we used in our simplified analysis. 
Nevertheless, we believe that our model provides a useful 
physical insight of the relevant aspects of the problem. 

In conclusion we have developed a powerful tool for 
the investigation of the interactions between condensates 
and we have demonstrated how they can quantitatively 
affect frequency, amplitude and shape of oscillations. In 
particular, the frequency shift measurement gives access 
to the number of atoms in the parent condensate N2 or, 
alternatively, to the (121 scattering length. This oppor- 
tunity seems particularly interesting for low N2, as we 
have found that the frequency shift scales roughly as 
log(iV 2 a 12 ) || . 

The agreement with the model, given its simplicity and 
lack of free parameters, is generally satisfactory. At least 



the discrepancy observed for the damping of center-of- 
mass oscillations suggests that the investigation of the 
relaxation beyond the mean-field theory (namely, by in- 
cluding the elastic scattering) could be an interesting sub- 
ject for future studies. The experimental perspectives of 
this system of two periodically overlapping condensates 
lead to the investigation of the interference between con- 
densates spatially separated during their evolution and, 
eventually, to studies of Josephson effects in two weakly 
linked condensates pl| j. 

We acknowledge fruitful discussions with F. Dalfovo. 
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FIG. 1. Center-of-mass oscillations as a function of the 
trapped evolution time (rf-pulse at t = 0) after 30 ms of 
ballistic expansion for (a) the non-interacting |1) condensate 
and (b) the two interacting condensates. The experimental 
points are compared with the relative theoretical predictions 
(discussed in the text). Each data point is the average of 
typically five measurements. 
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FIG. 2. Evolution of aspect ratio for the [2) condensate in 
the trap obtained solving the GP equations for the interact- 
ing (solid line) and non-interacting case (dashed line). In the 
interacting case the evolution of the aspect ratio is not a sim- 
ple superposition of the two oscillations with frequency uji ow 
and uihigh, as in the non-interacting case. We note also that 
the interaction between the two condensates produces an ini- 
tial delay for the onset of collective excitations, since the two 
condensates are created with the same density profile, and the 
scattering lengths 022 and 012 have the same value. 
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FIG. 3. Aspect ratio for the |2) condensate, as a function of 
the trapped evolution time, after 30 ms of ballistic expansion. 
The experimental points are compared with the theoretical 
predictions for the case of two interacting condensates. 
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